Quantum Dynamics of a Nanomagnet in a Rotating Field. 
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Quantum dynamics of a two-state spin system in a rotating magnetic field has been studied. 
Analytical and numerical results for the transition probability have been obtained along the lines of 
the Landau-Zener-Stueckelberg theory. The effect of various kinds of noise on the evolution of the 
system has been analyzed. 
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I. INTRODUCTION 



Molecular magnets with high spin and high mag- 
netic anisotropy can be prepared in long-living excited 
uantum states by simply applying a magnetic field 
lj. In a time dependent magnetic field, they exhibit 
stepwise magnetic hysteresis due to resonant quantum 
tunneling between spin levels Q. This phenomenon 
has been intensively studied theoretically within models 
employing Landau-Zener-Stueckelberg (LZS) effect 
til EL LI The formulation of theproblem, independently 
studied by Landau [3j, Zener and Stueckelberg |5| 
at the inception of quantum theory is this. Consider 
a system characterized by quantum states |1) and |2) 
with energies E\ and Ei respectively. Let the system be 
initially prepared in the lower-energy state, |1), and the 
field be changing such (due to, e.g., Zeeman interaction 
of the magnetic field with a spin) that E\ is shifting up 
while E2 is shifting down. After the levels cross and the 
distance between them continues to increase, the system, 
with LZS probability, P — exp(— ttA 2 /2hv), remains in 
the state |1). Here A is the tunnel splitting of |1) and |2) 
at the crossing, and v is the rate at which the energy bias 
between |1) and |2) is changing with time. This picture, 
of course, does not take into account any disturbance 
of the quantum states, |1) and |2), by the dissipative 
environment. Application of the conventional LZS effect 
_to_molecular magnets was initially suggested in Refs. 
Its dissipative counter par t was developed in Refs. 
13 IH 111 13 El 13 13 The amazing property 
of the LZS formula is that it is very robust against 
any effect of the environment 0, 0, 0|. This has 
allowed experimentalists to use the LZS expression to 
extract A in molecular magnets from bulk magnetization 
measurements El El HI El ■ 



The purpose of this paper is to elucidate the possibil- 
ity of a detailed study of quantum spin transitions and 
the effect of the environment in experiments with a ro- 
tating magnetic field. Quantum tunneling rates for a 
spin system in a rotating field have been studied before 
[22J. Here we are taking a different angle at this prob- 
lem, by computing the occupation numbers for quantum 
spin states. This approach can be useful for the descrip- 



tion of experiments that measure the time dependence 
of the magnetization. A weak high-frequency rotating 
field, Hq ~ 1 Oe, can be easily achieved electronically, 
by applying H x = H cos(wi) and H z = H sin(wt). The 
rotating field of large amplitude can be achieved by rotat- 
ing the sample in a constant magnetic field. In a typical 
molecular magnet, a rotating field not exceeding a few 
kOe will result in the crossing of two spin levels only, 
preserving the two-state approximation. We will com- 
pute the time evolution of the probability to occupy one 
of the two spin levels after a number of revolutions. We 
will demonstrate that this evolution depends crucially on 
whether the system is subject to the dissipative noise and 
that it depends strongly on the frequency of the noise. 
The equivalent of the LZS effect in the rotating field is 
studied in Sec. III. The effect of slow and fast noise on 
such probability in the rotating field is considered in Sec. 
IV . Consequences for experiment are discussed at the 
end of the paper. 



II. HAMILTONIAN 



We shall start with the Hamiltonian 
H = -DS 2 Z - g/j, B H ■ S . 



(1) 



where D is the uniaxial anisotropy constant, g is the 
gyromagnetic factor, fi B is the Bohr magneton, and H 
is the magnetic field, and S is an integer spin. We shall 
assume that the magnetic field rotates in the AZ-planc, 



H = H sm(ujt)e z + H cos(ujt)e x 



(2) 



At H = 0, the ground state of the model is double de- 
generate. The lowest energy states correspond to the 
parallel (m = S) and antiparallel (m = —S) orientation 
of S with respect to the anisotropy axis, with m being 
the magnetic quantum number for S. The effect of the 
external magnetic field is twofold. The Z-component of 
the field removes the degeneracy. The A-component pro- 
duces a term in the Hamiltonian that does not commute 
with S z . Consequently, at H x ^ the \m) states are no 
longer the eigenstates of the system. However, at 



gH B H<t: (25- 1)D 



(3) 



2 



we can treat the non-commuting term in the Hamiltonian 
as a perturbation. Throughout this article it will be as- 
sumed that the system is prepared initially in one of the 
saturated magnetic states, say |— 5) for certainty. This 
can be easily achieved at low temperature in molecular 
magnets with high easy-axis anisotropy. 

For the perturbation V(t) = gn B H ■ S the time- 
dependent perturbation theory gives the following ex- 
pression for the transition amplitude from the initial 
state, \—S), to any \m') state with m! ^ —S: 



Cm' (^) 



(m'\ gt jL B HS x \-S) x 



dtie 1 h cos(wii) 



J2 gn B HS* \m") (m"\g» B HS x |-5) 



dtie 1 * cos(ojti) x 



dt^e* h cos(wi 2 ) 



(4) 



where E m = —Dm 2 are the eigenstates of 7Yo = —DS 2 . 
When u> <C (E m > — E-s)/h~ for all ml , then the pertur- 
bation can be treated adiabatically. This requires the 
condition 
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FIG. 1: Time dependence of the energy levels (normalized 
by D) of the Hamiltonian (JSJ (solid line) at S = 4 and 
gfi B H/D — 3.1. Dash line shows the distance between the 
energy levels if they were unperturbed by the second term in 
Eq. ©. 



5, it decreases very fast as one moves away from the level 
crossing, that occurs at r„ = mr/2. This guarantees that 
the transitions between the two states are localized in 
time at the level crossing. This is clearly seen in Figure 
1 that shows the effect of the perturbation on the energy 
levels. 



hw < (25- 1)D, (5) 

that will be used throughout this paper. In molecular 
magnets, D is of the order of 1 K. Consequently, any to 
at or below GHz range satisfies Eq. JSJ. 

At T = 0, equations @ and (JSJ, and the initial con- 
dition, allow one to limit the consideration by the two 
lowest states, |5) and |— 5), of the unperturbed Hamil- 
tonian. Time-independent perturbation theory based 
upon the condition (0 permits the usual reduction of 
the spin Hamiltonian to the effective two-state Hamilto- 
nian for the tunnel split states originating from m = 5 
and m = —5, 

Heff = --h sm{r)a z + ^A(t)ct x . (6) 

Here a z , x are the Pauli matrices, h a = 1Sg[i B H is the 
amplitude of the energy bias, r = cot is dimensionless 
time, and A(r) is the tunnel splitting of |5) and |— 5) 
due to the transverse field H x {t) 0, 

A(r)=A cos 2S (r), (7) 

where 

_8S 2 D ( h \ 2S 

An important observation that follows from Eq. Q is 
that at T n = (2n + 1)tt/2 (with n = 0,±1,±2,...) the 
splitting is exactly zero. According to Eq. JJJ), at large 



III. DYNAMICS WITHOUT NOISE 
A. A single revolution 

In order to describe the evolution of the magnetiza- 
tion in a magnetic field rotating in the XZ-plane, we 
shall compute first the probability of staying at the ini- 
tial state after a rotation by 180 degrees, when a single 
level crossing takes place. 

The Schrodinger equation for the coefficients of the 
wave function, 

|*) =c_ s |-5)+c s |5) , (9) 
can be expressed as: 

„ d _ A(t) _ , , 

ihj t c s {t) = -h sin{ujt)cs(t) + ^c_ 5 (t) (10) 

where 5j(t) = Cj(t)exp[^ f t dt'h sin(u}t')]. In terms of 
the dimensionless variable, 

u =^—> (11) 
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Eq. lfTU|) turns into: 



. d 



cos (7u)cs(u) 



— sin(7w,)cs(u) + — cos 25 (7it)c_s (it) 
7 2 

(12) 



The problem is now denned by two dimensionless pa- 
rameters: 



A 2 
huj ho 



(13) 



which is similar to the parameter used in the LZS theory 
[IIH, and 



7 



Ao 
ho 



(14) 



which is a measure of the magnitude of the magnetic 
field. Notice that this parameter equals r c = u>t c , where 
t c = A /(uho) is the characteristic time of crossing the 
resonance. Thus, the condition 7 <C 1 is needed for the 
crossing to be well localized in time on the time scale of 
one revolution. This condition is required for the self- 
consistency of the method. 

From Eqs. 112|) one can compute numerically the time 
evolution of the coefficients c_s and cs, and, thus, the 
time evolution of the occupation numbers for any e, 7. In 
the fast (e <C 1) and slow (e ^> 1) rotation regimes, an- 
alytical formulas for the occupation probabilities can be 
obtained. These formulas are useful for further analysis. 
Inspection of Eas. (|12|) reveals that the deviation from the 
LZS result for slow rotation is small, since within the rel- 
evant time of the transition 5u ~ 1 and A(r) is nearly 
constant. On the contrary, for the fast rotation regime, 
(e <C 1), the relevant time interval is wider, Su ~ £ -1 / 2 . 
This allows a significant change of A during the tran- 
sition and makes possible a considerable deviation from 
the LZS result. 



1. Fast rotation (£<1) 

Following the procedure devised by Garanin and 
Schilling |2J], we can obtain the probability of staying 
at the initial state after a 180-degree rotation of the ex- 
ternal magnetic field. We choose the direction of the 
magnetic field to be initially antiparallel to the Z-axis. 
In the zero-th order of the perturbation theory c\ (u) = 1 . 
In the first order, such probability is then given by 



P = 1 



cos 2S (z) exp 



7 J- 



■ cos(z) 



dz 



(15) 



It is clear from this expression that for 7 < 1, only 
\z\ <C 1 contribute to the integral. This is in accordance 
with the fact that the transition takes place during the 
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FIG. 2: Probability of staying at |— S), P, as a function of 
parameter e at 7 = 0.1 and S = 5 . The solid line represents 
the numerical data, the short-dash line is the analytical re- 
sult, Eq. 1171 . and the long-dash line is the LZS result (see 
explanation in the text). 



time interval Ao/houj, which is of order 7 compared with 
the time of the integration. Consequently, one can ap- 
proximate cos(z) by 



cos(z) 



1 



2 ' 



cos\z) 



e 22 , 



(16) 



where the exponential form is chosen to insure fast con- 
vergence of the integral, and set infinite integration limits 
in Eq. (JTSJ: 



P 



e 

7 J-c 



7T 

l -2 



exp 

> 



+ S)z 



dz 



(17) 



The probability, P, of staying at the initial state \—S) 
after a 180-degree rotation is shown in Fig. El for differ- 
ent values of e. In the figure, the numerical results are 
compared with the result given by Eq. (|T7|) , and with the 
LZS result, 



Plzs = exp[-7re/2] . 



(18) 



Note that Eq. (|18|l would be our result for the probability 
if we applied the LZS theory to the version of the Hamil- 
tonian © that is linearized on t. As can be seen from 
Fig. at 7, e < 1, Eq. (|T7|) provides a good approxima- 
tion. The difference between the numerical result and the 
LZS result is considerable. Fig.|^lshows P for different 7 
at e = 0.1. Here again Eq. I|17fl . but not the LZS formula, 
provides a good approximation for the 7-dependence of 
the staying probability. 



2. Slow rotation (e>lj 

In the case of a slow rotation it is convenient to seek the 
solution of the Schrodinger's equation in the adiabatic 
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FIG. 3: Probability of staying at the initial state, P, as a 
function of 7 at e = 0.1 and S = 5 . The solid line represents 
numerical data, the short-dash line is the analytical result 
Eq. 11711 . and the long-dash line is the LZS result. 



FIG. 4: Probability of staying at the initial state, P, as a 
function of the parameter e at 7 = 0.1 and S — 5 . The solid 
line represents the numerical data, the short-dash line is the 
analytical result, Eq. 1271 . and the long-dash line is the LZS 
result 



basis of the two-state Hamiltonian. Then, one can follow 
a procedure similar to that used for the fast rotation 
regime. The adiabatic basis is given by: 



\*±) = -^(±k ± \-S) + k T \S)) 



where 



fct(r) = Jl± 



W{t) 



^1F 2 (t) + A 2 (t) 



(19) 



(20) 



with W(t) = /i sin(r) and A(r) = A cos 26 (t). The 
corresponding adiabatic energy levels are 



1 



E±(t) = ±-v/^ 2 (r) + A2(r). 

Expressing the wave function as 

|*)=c+|* + )+c_|*_) 
we can now write the Schrodinger equation for 



c±(t) = exp 



i J dtE_(t) 



c±(t) 



(21) 



(22) 



(23) 



in terms of the dimensionless variable u defined above: 
d . 



du 
d_ 
du 



N _ 1 dw/du _ 
— ieil('ju)c + — — — — C- 



2 (1 + w 2 ) 



c_ = — 



1 dw I du _ 
2(l + w 2 ) C+ 



where 



— sin 2 (z) + cos 4S (z) 



(24) 



(25) 



and w(z) = W(z)/A(z). 



The asymptotic behavior of c_ (u) is c_ (u) — > 1 for 
e —> 00. At e 3> 1 the coefficient c_(u) remains close to 
1. Also, far from the crossing point, c + (u) w c_5, as can 
be seen from Eqs. I|19|l and (|2()|l . These two facts allow 
one to obtain the probability of staying at the initial state 
in the first order of perturbation theory: 



P = 



1 



dz 



dw/dz 
(l + w 2 ) 



exp 



i- 1 dz'n(z') 
7 Jo 



(26) 



At 7 1 the integral is dominated by z close to zero. 
The correct prefactor can be obtained by applying the 
procedure outlined in Ref. |3 With the accuracy to 7 2 
this gives: 



P 



exp 



tic 

Y 



1 



5'- 



1 



7 



(27) 



This probability is shown in Fig. 0] for different values of 
e. The numerical results are compared in the figure with 
the result given by Eq. (|27|l . and with the LZS result. 
As can be seen from Fig. at 7 <C 1, Eq. (|2*7|l provides 
a good approximation for e > 5. 



B. Continuous Rotation 

Our treatment of the continuous rotation is based upon 
the smallness of the transition time in comparison with 
the period of the rotation. Periodically driven two-state 
systems of that kind have been studied before |14| . 

The individual crossing is described by the transfer 
matrix: 



P 



M = 



(28) 



5 



where 9 is the Stokes phase given by: 



4 + arg 



t r H): 


e 






+ 4 





(29) 



This matrix transforms a given initial state into the 
after-crossing final state in terms of the unperturbed ba- 
sis. The above expression corresponds to the crossing 
in which the \—S) level moves up towards the \S) level 
that is moving down. In the opposite case, M should be 
replaced by the transpose matrix, M T . 

To describe the evolution of the system between cross- 
ings it must be noted that, as we have shown before, 
far from the crossings the unperturbed basis \—S) , \S) 
almost coincides with the adiabatic basis |+),|— ), see 
equations (|T§jl and lj2"U|l . In this region, the evolutions of 
|— S) and | S) are then considered independent, so that 
they can be described by the propagator 



exp[(- 



ii+i 



ia] 



exp[(— l) n ia] 



(30) 




FIG. 5: Time dependence of the probability of finding a 
continuously rotating system in the initial state \—S) for 
e — 0.04, 7 = 0.06 and S = 5. The solid line represents nu- 
merical data. The dash line shows analytical result obtained 
by successive application of M, M T and G n ■ 



where 



1 / D 



*=- 2 {j-)J o drVWWTAW. (31) 

With the help of Eq. (|2"51) and Eq. l(3TJ|) one can compute 
the time evolution of the coefficients c_s,cs in Eq. 0. 
Starting with the initial state, the state of the system 
after the n-th crossing can be obtained by the succes- 
sive action of M,M T and G n . The time-dependence of 
the probability of finding a continuously rotating system 
in the initial state |— S) is shown in Fig. [SJ The figure 
shows good agreement of the above analytical method 
with numerical calculation. It is important to notice that 
in the absence of dissipation the system does not arrive 
to any asymptotic state at r — > oo. The behavior of 
the probability shows a long-term memory of the initial 
state, which is somewhat surprising. This prediction of 
the theory can be tested in real experiment. 
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FIG. 6: Time dependence of the probability of finding a par- 
ticle in the state \ —S) in the presence of a low- frequency noise 
for e = 1/3,5" = 10, and 7 = 0.01. The plotted probability 
is the average over an ensemble of 20 two-state dissipative 
systems. 



IV. DYNAMICS WITH NOISE. 

When considering the effect of the noise, it is important 
to distinguish between the following three regimes: 

7 < 1< lj/T (32) 
7 < w/r < 1 (33) 
lo/T < 7 < 1 (34) 

where T is the characteristic frequency of the noise. The 
first of these conditions corresponds to the situation when 
a few revolutions may occur before any contribution of 
the noise becomes apparent. Consequently, during the 
time interval satisfying t < 1/T one can use the results 
for the probability obtained in the previous section. Un- 
der the condition (|33|l . one can use the previously ob- 
tained results for a singular crossing but needs to take 



into account the destruction of the relative phase of the 
two states by the noise before the next crossing takes 
place. Under the condition <|51|l the results of the pre- 
vious section do not apply because the coherence of the 
quantum state is destroyed by the noise on a timescale 
that is less than the crossing time. 



A. Low- frequency noise, 7 <C lo/T <§C 1 

The situation corresponding to the condition (|33|) can 
be easily described by the Hamiltonian 

Heff = ~h sin(r)cr z + ~A(t)ct x - i]<r z , (35) 



6 



where 7/(t) <C ho is a random magnetic field in the in- 
direction, with the correlator 



(r,(r)r 1 (T'))=rile[(uj/V) 



(36) 



© being the theta-function. 

The time dependence of the probability of finding 
the system in the state |— 5) is shown in Fig. © for 
e = 1/3, 5 = 10, and 7 = 0.01. In accordance with the 
expectation, the probability for the system to occupy the 
state with m = —S (or m = 5), after going through a 
few oscillations, tends to the asymptotic value of 0.5. For 
molecular magnets, only the average of the probability 
over an ensemble of two-state dissipative systems is of 
practical importance. The probability shown in Fig. © 
is such an average. 
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FIG. 7: Time dependence of the probability of finding a 
particle in the state |— S) in the presence of a high-frequency 
noise for S = 10, e = 1/3 and i = 0.01 . 



High-frequency noise, uj/T <C 7 <C 1 



In this limit the coherence is completely suppressed, 
by, e.g., interaction with phonons, and the evolution of 
the population of energy levels must be described by the 
density matrix. In this case the population iV_s, of the 
initially occupied state |— 5) is given by (l3l.l25| 



s 



AW 2 

dt 2 

The solution is 



hT 



w{t) 2 + (hry 



;(N_s~N s ). (37) 



'2 9 



N-s = ~ jl + exp 

where g(z; b; a) is given by 
g(z; b; a) — 

(2n + 2)/ 1 (|)-/ 1 (|-f 



r;25;^ 
h 



(38) 



(39) 

if z = (2n+l)| + £ 



g(z;b;a)=2nf 1 (-) + f 1 (0, ifz = n7r + £, (40) 



sin 2 (z)\ sin(z) 



where n = 0, 1, 2..., < £ < tt/2, 
/iW = fi Q;i-M;|sin 2 (z), 

(41) 

and Fi(a,bi,b2',c;x,y) is the Appell hypergeometric 
function of two variables. 

The time-dependence of the occupation of the state 
is shown in Fig. dfor e = 1/3,5= 10, hT/h = 0.01. 
As in the case of a low-frequency noise, the probability 
to occupy either of the two levels tends to 0.5 after a 
few revolutions. The difference between the two regimes 
is that in the case of the high-frequency noise the prob- 
ability monotonically approaches the asymptotic value 
without exhibiting any oscillation. 



V. CONCLUSIONS 



We have studied the equivalent of the LZS effect for a 
spin system in a rotating magnetic field. Typical time de- 
pendence of the probability of staying at the initial state 
has been computed for three different situations. The 
first is the situation when the noise is irrelevant on the 
time scale of the measurement, Fig. |SJ In this case the 
system exhibits coherent behavior and long-term mem- 
ory effects. The second situation corresponds to the noise 
that decohere quantum states within the time of each 
revolution but is slow enough to provide pure quantum 
dynamics during the level crossing, Fig. The third 
situation corresponds to a very fast noise that does not 
allow the use of wave functions for the description of the 
crossing and requires the density-matrix formalism, Fig. 
When the noise becomes important, the occupation 
probability of each level approaches 1/2 after several rev- 
olutions. However, the asymptotic behavior depends on 
the frequency of the noise. The three regimes discussed 
above are given by equations l|3^|) . l|3*3*|l . (|3"4l . One must be 
able to switch between different regimes by changing the 
angular velocity of the rotating field and/or temperature. 
Experiments of that kind can shed light on the effect of 
dissipative environment on the resonant spin tunneling in 
molecular magnets. To be on a cautious side, one should 
notice that the evolution of the magnetization in a crystal 
of magnetic molecules also depends on the dipolar inter- 
actions between the molecules [2(| H3, • Our results 
are likely to be relevant to molecular magnets when the 
amplitude of the rotating field significantly exceeds dipo- 
lar fields. A candidate for such a study would be, e.g., 
an uniaxial Ni-4 molecular magnet that has no nuclear 
spins and that, with good accuracy, is described by the 
Hamiltonian studied in this paper. 
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